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Further investigation of the participant plane correlations within a Glauber model framework is 
presented, focusing on correlations between three or four participant planes of different order. A 
strong correlation is observed for cos(2<£>2 + 3<&3 — 5$s) which is a reflection of the elliptic shape of the 
overlap region. The correlation between the corresponding experimental reaction plane angles can 
be easily measured. Strong correlations of similar geometric origin are also observed for cos(2$2 + 
4$J - 6$g), cos(2$2 - 3$3 - 4<E>J + 5<E>s), cos(6$2 + 3$3 - 4$J - 53>s), cos($J - 2$2 - 3$3 + 4$J), 
cos($i +6$2 — 3<&3 — 4$J), and cos($J; +2<E>2 +3$3 — 6^e), which are also measurable. Experimental 
measurements of the corresponding reaction plane correlators in heavy ion collisions at RHIC and 
the LHC may improve our understanding of the physics underlying the measured higher order flow 
harmonics. 
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In a previous paper [1], one of us proposed a method 
for measuring the correlations between several reaction 
planes of different order. Wc estimated the magnitude 
of these correlations in configuration space via a Monte 
Carlo Glauber model, and several strong spatial correla- 
tors were identified. Despite the possible non-linear mix- 
ing between harmonics of different order in the hydrody- 
namic evolution [2-4], these geometric correlations may 
still survive and contribute to the reaction plane correla- 
tions in momentum space. In this paper we discuss sev- 
eral geometric correlators involving three and four par- 
ticipant planes, which are of current experimental inter- 
est, and which are not covered in the literature. Related 
correlators have been studied both numerically [5-7] and 
analytically [9]. 

As pointed out in Ref. [7, 8], the reaction plane cor- 
relations that can be measured experimentally involve 
various linear combinations of the n-th order planes $„, 
ci$i + 2c 2 < E>2-- - + lci$i, where the integers (ci . . . cj) sat- 
isfy the constraint 



ci + 2c : 



2 — 



lci = 



(1) 



Due to this constraint, only I — 1 angles are independent. 
The differential distribution in this observable is an even 
function and can be expanded into a Fourier series 



d(ci$i + ... + lci®i) 



3=1 

COSj(ci$i + ... +lci$i) , 



where 
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The Fourier coefficients can be determined from the ex- 
perimentally measured event plane angle ^ n and associ- 
ated resolution factor Res{jc n n,^I'„} 
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Res{jci*i}...Res{jc^*J 
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where 



Res{jc n nV n } = (cos jc n n(^ n - $„)) . 



(5) 



The precision with which these reaction plane correla- 
tions can be measured is limited by the magnitude of the 
resolution, which is expressed in terms of the resolution 
parameter \ n [10]: 



Reslmn^n} 
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I ( m -l)/2(-^) + I( m +l)/2(-^-) 
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In general, Xn and hence Rcs{mn^! n } decrease quickly 
for increasing n. For event plane measured in 3 < \r)\ < 5 
in Pb+Pb collisions at the LHC, the ATLAS Collabora- 
tion shows that Xn decreases from about 2 for n=2 to 
about 0.08 for n = 6, and is negligible for n > 7 [11]. In 
contrast, ~R.es {runoff „} decreases more slowly with m at 
fixed n, especially for n=2 and 3 cases where Xn is close 
to unity. The dependence of H,es{mn^ n } on n and m 
limits the types of correlations that are accessible to the 
experiments. 

The three-plane correlator can be generally expressed 
as a linear combination of two two-plane correlators 



c„n<I>„ + c m m$ ri 



cil$i = c m m($ m - $„) 

+ al(<S>i - $„) 



which is redefined in terms of $ a .b = ($ a — 

c„n$„ + c m m$ m + Cll®l = C m ®m,n + Cl$l, 



(7) 
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Here n < m < I, and we have used the constraint in 
Eq. 1. We shall refer to these three plane correlations as 
"1-m-n" correlations. 

The correlation signals can be accessed via a Fourier 
expansion of the event distribution in 1 (& m ,n, $i,n)'- 



oc 



1 + 2 V i,n COS j$ m ,n + Vj> n COS J$, )B 



+ 2 E ^m,! C0S (**m,» ± ■ (8) 

The meaningful coefficients arc those that satisfy the con- 
straint of Eq. 1: (jm mod n) = for V^, „, (j7 mod n) 

for Vjj n , and (im ± jl mod ra) = for V/'^. 

The discussion so far involves experimentally measured 
correlations between reaction plane angles $„ or event 
plane angles ty n , which are defined in the momentum 
space. These correlations are partially related to analo- 
gous correlations between participant planes in the 
initial geometry. Previous study of participant plane 
correlations focused on three-plane correlators contain- 
ing $|. In this work, we explore three-plane correlators 
that do not involve as well as various four-plane cor- 
relators. These correlations are estimated with Monte 
Carlo Glauber simulations of Au+Au collisions using a 
nucleon-nucleon cross-section of a = 42 mb [12]. The 
and the eccentricity, e„ , are defined through the distribu- 
tion of participants and binary collisions in the transverse 
plane, with a weight of 8 = 0.14 for binary collisions, and 
(1 -S)/2 = 0.43 for participants [14] 



.» _ (r 2 e in ' 1 ') 
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Here (r, <j)) are measured relative to the weighted center 
of mass [13]. Alternatively, 3>* can be defined with an 
r 3 -weight for n = 1, and a r n -wcight for n > 1, and this 
definition is referred to as r n -weighting. We also calcu- 
lated the participant plane angles with CGC simulations 
using both the r 2 and the r™-weighting [15]. Finally, we 
note that <I>* defines the major axes of the eccentricity [1] 
and is rotated by ir/n relative to traditional definition 
based on the minor axis. 

Two interesting three-plane correlators are, 



c 2 $2 + 3c 3 $J + 5c 5 $; = c 3 <f>l 2 + c 5 $ 



-5,2 i 



(10) 



and 



c 2 $; + 4c 4 $4 + 6c 6 $g = c 4 $4 :2 + c 6 $* 2 . (11) 



1 This expression can be obtained from a double Fourier series in- 
volving cosi<J> mi „ cos j<E>( n and sini<3? m>n sinj&i n and, in prin- 
ciple, the corresponding mixed terms. However, terms lin- 
ear in sine vanish since the event distribution is even under, 

*m.n -> -<J>m.n. 



Fig. 1 summarizes the "5-3-2" correlations present in 
Glauber and CGC simulations for the two different 
weighting schemes. A strong diagonal correlation is ob- 
served, corresponding to = or (cosi(2$2 + 
34>3 — 5*l>5)) . The coefficients are nearly zero for other val- 
ues of (i, j). The origin of this correlation is similar to the 
well known "3-2-1" correlation, (cos($|+2$?;-3$3)) [6]. 

Specifically, both of these correlations are geometric 
and of order t2- To see this we will use a simplified (but 
less accurate) version of the independent cluster model 
[9] and compute a closely related correlation function 



(6 3 e 5 cos(2$; + 3$g - 5$g)) 



(12) 



In the cluster model, N independent clusters are drawn 
from a distribution n(x), which is the average number of 
clusters per unit area in the transverse plane, with x = 
(x,y). n(x) is proportional to the participant density 
in an optical Glauber model. The cluster density in a 
specific event is n{x) = 5n(x)+n(x) where 8n(x) satisfies 
Poisson statistics 



(5n(x)5n(y)) = n(x)S 2 (x — y) 



(13) 



For a given event, the odd participant planes and their 
angles are given by the integral 



m(tf>* — <b* R ) 



1 



N{r 2 ) 



^xSnix)^ 1 ^-^ , 



(14) 

where the r x and <p x are the radius and azimuthal angle of 
the cluster, and we are measuring all angles with respect 
to the original reaction plane $>* R ~ $2 + 7r /2 ) and are 
working to leading order in 1/N. Multiplying Eq. 14 by 
its conjugate and averaging over the statistics of 5n[x), 
we find 



(4) 



N(r 2 Y 



(15) 



which explains (again) why e n does not decrease with 
n in the Glauber model [9]. Similarly, after construct- 
ing e3£5e^ 2 *« +3 * 3_5 * 5 ^ with Eq. 14, and averaging over 
Sn(x), we find 

(6 5 e 3 cos(2$; + 3$g - 5$g)) (r 4 cos(2(^ x - 



(r A ) 



(16) 



Eq. 16 clearly shows the geometric origin of these corre- 
lations. In fact, all "(n+2)-n-2" correlations with n odd 
are equal, and given by Eq. 16. When the independent 
source model is improved by including the shift in the 
center of mass and terms suppressed by 1/N, it quanti- 
tatively describes the results of full Monte-Carlo Glauber 
and CGC simulations [9]. 

A similar correlation in the Glauber model is studied in 
Fig. 2 which examines the "6-4-2" three plane correlators. 
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FIG. 1. (Color online) The normalized event distribution for 3($3 — $2) and 5(<E>5 — $2), i.e. d 2 A r cvts /[d$| 2 d^s^] (top row). 
The corresponding Fourier coefficients (bottom row) in 40-50% centrality class, see Eq. 8. Note that the constraint in Eq. 1 
requires i,j to be both even or both odd. 
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FIG. 2. (Color online) The normalized event distribution for 4(<E>4 — $2) and 6($6 — $2)1 *- e - d 2 A jr e vts/[d$2 :2 d$6,2] (top row) 
The corresponding Fourier coefficients (bottom row) in 40-50% centrality class, see Eq. 8. 



The largest term corresponding to (cos(2$2+4<f>|-6<f>g)), 
is much bigger than (cos6($g — $2)) an d almost as big 
as (cos4($| — $2))- This term captures the strong pos- 
itive correlation of 4>g and $| relative to the $2 plane 
as shown in the top panels of Fig. 2. However, in con- 
trast to the "5-3-2" correlation in Fig. 1, the alignment 



of <£>g, 3>| relative to $2 i g not strictly along the di- 
agonal ($g j2 ~ 1.7$| 2 ), leading to large coefficients 
along = + 1) and (i, — i + 2). This 

non-diagonal behavior reflects a strong influence of aver- 
age geometry to even-order participant planes $|,$g..., 
which leads to strong correlations between $2n an£ l 





FIG. 3. (Color online) The centrality dependence of several 
three-plane correlators. 



FIG. 4. (Color online) The centrality dependence of several 
four-plane correlators not involving <E>i plane. 



and hence narrow peak in distributions of $| 2 an d 
$g 2 [1]. After removing these two-plane correlation com- 
ponents 



d 2 N r 



dAUts dN r 



m, 2 



d$ 



(17) 
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the resulting 2-D distribution is modulated around di- 
agonal direction and with a similar magnitude as in the 
"5-3-2" case. This harmonic variation around the diag- 
onal band (also see Fig 1(b) and (d)) can be analyzed 
similar to the "1-2-3" case [6]. 

Figure 3 shows the centrality dependence of several 
three-plane correlators for which the corresponding ex- 
perimental event planes may still have decent resolution. 
Very strong signals are observed for (cos(2$2 + 33>jj — 
5$s)) and (cos(2$2 + 4$J - 6$g)); the signals are even 
bigger than (cos($J + 2$2 — 3$g)) and are comparable 
to (cos4(^-$|)). The value of (cos(2$|-8$| + 6$|)) 
is also large. In contrast, the values for other correlators 
are small or even slightly negative for the r 2 -weighting in 
mid-central collisions. 

The behavior of these correlators towards central col- 
lisions is also quite interesting. This is the region where 
all the correlations are presumably dominated by fluctua- 
tions of participating nucleons. The values of (cos(2$2 + 
3$^ - and (cos(2$^ +4$| - 6$£)) seem to reach a 

constant at 10% for iV par t > 300, at least for the Glauber 
model. Similar observations are also made previously [1] 
for (cos4($ 2 ! - $2)), (cos6($^ - <&%)), (cos2($? - 
(cos($ 1 +2$^-3$^)) and (cos($£ +3$£ - 4$|)). Since 
the non-linear mixing between different harmonics due 



to hydrodynamic evolution is expected to be relatively 
small in central collisions, measuring the corresponding 
reaction plane correlators in this region may provide some 
handle on the relative role of the linear and non- linear re- 
sponse [3, 4] 

The four-plane correlator in the Glauber and CGC 
models can be analyzed with a Fourier analysis similar 
to Eq. 8. Instead of presenting a general analysis, we 
will simply discuss several participant plane correlators 
that have a large signal and acceptable reaction plane 
resolution for the corresponding momentum space mea- 
surements. 

Figure 4 shows the several four-plane correlators which 
do not involve the §\ plane. The strong signal observed 
for four of the correlators shown can be understood as the 
inter-correlation between two correlators involving two or 
three planes, which each have strong signal. For example, 
the first pair of four-plane correlators 



2$; - 3$; - 4$J + 5$5 =4($a - K) - (2$2 + 3$a 
6$2 + 3"1>3 - 4$2 - 5$5 =4($2 - $4) + (2$| + 3$3 



5$s) , 
(18) 



are very strong (the black circle and blue cross symbols), 
while second pair 



$2 + 6$; - 4$: - 6$; =4($* - $t) - 6($* - $*) 

6(*5-$$) : 



4$; - 6$; - 4$2 + 6$; =4($^ - <f>* 4 



(19) 



are somewhat smaller (the green diamond and magenta 
square symbols) since the 6(<E>3 — <E>g) correlator does not 
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FIG. 5. (Color online) The centrality dependence of several 
four-plane correlators involving *i plane. 



source the elliptic shape. Furthermore, since 



2(sin4(* 2 ' - *J) sin(2«I>^ + 3*3 - 5*£)) = 
(cos(2*;!; - 3*3 - 4*| + 5*5)) - (cos(6*2 + 3*3 

2{sin4(<I>2 - *1) sin6($3 - *g)) = 
(cos(*2 + - 4*J - 6*6)) - (cos(4*2 - 6*3 ■ 



4*2 



1*1 



-5*5)> 

-6*|)> 
(20) 



each pair of correlators also allow us to infer the relative 
sign of the two composing correlators. For example, the 
centrality dependence of the relative magnitude of both 
pair of correlators in Fig. 4 suggest that the average val- 
ues of these sine products are positive in mid-central and 



peripheral collisions for r" -weighting, while they remain 
negative for r 2 -weighting. The last four-plane correlator 
(the red star symbols) reflects inter-correlation between 
two three-plane correlators 



4$* 



3^ 



(2<5>* 2 - 6$; + 4$*) + (2$; + 3$3 - 5$^) , (21) 



and is small. 

Lastly, Fig. 5 shows five four-plane correlators that 
have a large signal and contain the §1 plane. These 
curves reflect the correlation between a three-plane cor- 
relator and a two-plane correlator 



$1 



2$5 - 3$3 + 4$4 = ($t + 2$2 
6$| - - 4$2 = ($1 + 2& 2 
2<E>! + 3$3 - 6$6 = ($* + 2$2 



3*3) 
3$3) 

ml) 



4($ 2 



$* _ 2$* _ 4$2 + 5$5 = ($* - 6<E>2 + 5<E>5) - 
3$* - 2$a + 3$3 - 4$2 = -($* + 2$| - 3$3 



$1) 

*S) 
*!) 

- 4($2 - $2) 

+ 4($* - $2) 
(22) 

Since the two composing correlators each have strong 
signals and are correlated with either the $2 -plane (the 
first four) or the plane (the last one), it is not sur- 
prising that these correlators also have sizable signal in 
mid-central collisions. 

In summary, correlations involving three or four partic- 
ipant planes are investigated in a Glauber model frame- 
work. These correlations are calculated in the configura- 
tion space, but are expected to contribute to the event 
plane correlations in momentum space, especially in cen- 
tral collisions. Several significant correlators are identi- 
fied and the reason for their large magnitudes are clari- 
fied. Many of these correlators are expected to have de- 
cent resolutions in Au+Au or Pb+Pb collisions at RHIC 
and the LHC, so should be measurable if the signal are 
as big as predicted by the Glauber model. 
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